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Abstract 



\ We demonstrate that for a finite-size quark-gluon plasma the induced gluon radiation 
from heavy quarks is stronger than that for light quarks when the gluon formation length 
becomes comparable with (or exceeds) the size of the plasma. The effect is due to oscilla- 
tions of the light-cone wave function for the in-medium q — ► gq transition. The dead cone 
Q-f model by Dokshitzer and Kharzeev neglecting quantum finite-size effects is not valid in 
this regime. The finite-size effects also enhance the photon emission from heavy quarks. 

1. It is widely believed that the observed at RHIC strong suppression of high-p^ hadrons 
(jet quenching) (for a review, see e.g. [1]) is due to radiative and collisional parton energy 

j> ■ loss in the quark-gluon plasma (QGP) produced in the initial stage of AA-collision. The 
dominating contribution to the energy loss comes from the induced gluon emission (for 
reviews, see e.g. [2, 3]). The effect of collisional loss is relatively small [4, 5, 6]. One of 
the interesting questions, which is important for jet tomography of /1/1-collisions at RHIC 
and LHC, is the question on the difference between gluon emission from light and heavy 
quarks. Dokshitzer and Kharzeev [7] suggested that gluon radiation from heavy quarks 
should be suppressed due to the dead cone effect. However, the RHIC data on suppression 

J> ! of the non-photonic electrons [8] indicate that the energy loss of heavy quarks may be 
similar to that for light quarks. The theoretical calculations of the energy loss of one of 
us [5] within the light-cone path integral (LCPI) approach [9, 10] also are in contradiction 
with the predictions of [7] since for a finite-size (FS) QGP at high energies the radiative 
loss has an anomalous mass dependence, i.e., AEheavy > A£^ 9 ta- 

In the present paper we give a physical interpretation of the anomalous mass depen- 
dence of the induced gluon emission. We show that the effect is related to the quantum 
FS effects which come into play when the gluon formation time becomes comparable with 
the size of the QGP. Physically the anomalous mass dependence of the induced gluon 
emission is due to oscillations of the light-cone wave function (LCWF) for the in-medium 
q — > gq transition. The present paper restricts its attention to the physical nature of the 
effect. For this reason, to make the analysis more transparent we, similarly to [7], consider 
a QGP with a constant density. The mass dependence of the jet quenching parameter for 
expanding QGP will be addressed elsewhere. 

2. We consider a fast quark with energy E produced in a QGP at z = (we choose the 
z axis along the quark momentum). In the LCPI formalism [9, 10], which we use, the 
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probability of the induced gluon emission can be written in terms of the Green's function 
for a two-dimensional Schrodinger equation with the Hamiltonian 

1 / <9\ 2 in(z)a 3 (p,x) 1 

H = -^[o- p ) 2 + v (1) 

where x is the gluon fractional momentum, fx = Ex(l — x), n is the number medium 
density, 0-3 is the cross section of interaction with a plasma constituent of the qqg system, 
Lf = 2/i/e 2 , e 2 = m 2 q x 2 + m 2 g {l — x) (m q and m g are the quark and gluon quasiparti- 
cle masses). In the low density limit Lf gives the coherence/formation length of gluon 
emission in an infinite medium (the Bethe-Heitler regime). The three-body cross section 
entering (1) can be written as a 3 (p,x) = {9[<7 2 (p) + cr 2 ((l ~~ x )p)] ~ °"2(^p)}/8, where 
a 2 (p) is the dipole cross section of interaction with color center of the color singlet qq 
system. The Hamiltonian (1) describes in- medium evolution of the LCWF of a fictitious 
qqg system. 

The gluon spectrum for a FS medium can also be written in the form [11, 12] 



-dx- = ! dZn{z) dx—' (2) 



do {x z) f 

eff dx ! — = Re J dp^*(p,x)a 3 (p,x)tp(p,x,z) , (3) 

where L is the quark pathlength in the medium, ip(p, x) is the ordinary LCWF for q — > gq 
transition in vacuum, and ip(p,x,z) is the in-medium LCWF at the longitudinal coordi- 
nate z (hereafter we drop spin and color indices). The latter can be written as 



z 

ip(x, p,z) = U J dzlC(p, z\p',0) 



(4) 

p'=o 



where K, is the Green's function for the Hamiltonian (1), and U is the spin vertex operator 
(its specific form can be found e.g. in [10, 13]). In the low density limit tp(p, x, z) goes to 
ip(p,x) at z/ Lf ^> 1. In this limit (3) reduces to the ordinary Bethe-Heitler cross section, 
dcr BH /dx, in the form derived in [14]. 

3. Let us consider the qualitative pattern of gluon emission. The integral (4) is saturated 
at z J> Lf, where Lf is the effective (with the LPM effect) gluon formation length in 
an infinite medium. The typical transverse size of the qqg system for gluon emission in 
this regime, p, is related to Lf by the Schrodinger diffusion relation p 2 ~ Lf/ fx. In a FS 
medium the dynamics of the gluon emission depends crucially on the ratio £ = L/Lf. 
For £ ^> 1 ip(p,x,z) is very close to ip(p, x, 00). In this regime the FS effects become 
negligible and the spectrum is close to that for an infinite medium (we call this situation 
the infinite medium regime). At £ < 1 the effective Bethe-Heitler cross section is chiefly 
controlled by the FS effects. In this regime (as in [12] we call it the diffusion regime) 
the dominating contribution comes from N — 1 rescattering [12, 15, 16] which gives the 
effective Bethe-Heitler cross section of the order of ^da BH /dx. 
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In the infinite medium regime the Coulomb effects are not very important and the 
gluon yield can be estimated in the oscillator approximation (OA) which corresponds to 
the parametrization cr 2 (p) = C 2 p 2 . Then the Hamiltonian takes the oscillator form with 
the oscillator frequency Q = \/—inC 3 / \i with C3 = | {9[1 + (1 — x) 2 } — x 2 } C 2 • Note that 
in terms of the well known BDMPS transport coefficient q [17, 2] C 2 = gCV/^nC^. In 
the OA the probability of gluon emission per unit length reads 



dP 
dxdL 



drr BH 



(5) 



where da^/dx = 2a s C 3 P gq {x) fine 2 , is the Bethe-Heitler cross section (P gq (x) is the 
ordinary q — > g splitting function), and S{rj) is the LPM suppression factor given by [10] 



S(V) 



sh y 



cxp 



VV2, 



cos 



+ sin 



(6) 



with r] = Lf\Q\ = ^AnC 3 p/e 2 . At r] 3> 1 when the LPM suppression becomes strong 
from (6) one can obtain S(rj) ~ {\ — -^75)- Then from (5) one obtains 



dP 
dxdL 



2^ 



j2r,C 3 


[l- ^ 1 







(7) 



Using (7) we obtain the heavy-to-light mass suppression if-factor 

n(M 2 - m 2 ^ 2 



K 



Ay/E(l - x)nC 3 ' 



(8) 



where Mq is the heavy quark mass. Thus, we see that, similarly to the ordinary Bethe- 
Heitler cross section, in the infinite medium regime with strong LPM suppression the 
gluon yield falls with quark mass. Note that the analysis [7] is supposed to be applicable 
namely to the infinite medium regime. In our notation the mass fT-factor obtained in [7] 
reads (in [7] the light quark and gluon are massless) 



K DK = 



1 + 



M 2 x 3 / 2 
3^/EnC 2 /2 



(9) 



Both (8) and (9) are obtained for strong LPM suppression using the OA. Some difference 
between the two expressions is not surprising since our formula, contrary to that of [7], 
is obtained with accurate treatment of the mass dependence of the oscillator Green's 
function. 

Let us now consider the diffusion situation. For a given L the boundary x for onset of 
this regime can be obtained using the OA estimate Lf ~ min(Lj, [fi^ 1 ) [12]. Assuming 
that the LPM effect is strong in the infinite medium regime (it means that Lf ~ |^| _1 ) 
one obtains that the diffusion situation corresponds to x ^nC^L 2 / E (we assume that x 
is small). Of course, it is only a crude estimate, and as will be seen from our numerical 
results in reality the FS effects becomes important at much smaller x. 
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Neglecting a small contribution from N > 2 rescatterings the effective cross section 
(3) for Debye potential with screening mass m D can be written in momentum space as 
[12, 18] 

da^ f H f {x,z) a 3 s C T P gq (x) \ , 771/-1 \ woi 

[F(l, z) + F(l — x, z) — F(x, z) /9J , 



n 2 C F 
dkdp 

(k 2 + m 2 D ) 2 ' 



H(yk,p)-[l-cos((p 2 + e 2 )p 2 d (z)) 



H(k,p) 



(P - k)p 



(p 2 + e 2 ) 2 (p 2 + e 2 )((p-k) 2 + e 2 ) ' 



(10) 

(11) 
(12) 



where Ct,f are the plasma constituent and quark Casimirs, Pd(z) = y z/2/j, is the diffusion 
radius. We represent F as a sum F = F + SF, where F Q = F(e = 0), and 5F is the 
mass correction. In the massless limit the momentum integration in (11) gives F (y, z) = 
7T 3 y 2 p 2 ,(z)/2. This leads to the spectrum [12] 



dP, 



N=l 



dx 



imL 2 aiC T P 9q (x)\l + (1 



X) 



x 2 /9] 



e=0 



8C F Ex(l - x) 



(13) 



An exact analytical calculation of 5F is impossible. We have performed approximate 
calculation of 5F for £ ^> 1. Keeping only the terms with large logarithms one can obtain 
in this limit 



my, ,) » U >og 2 ( ^rr) + log ( ') log / ' 2 



(14) 



Then, neglecting in (14) the linear subleading logarithms, we obtain for the mass correc- 
tion to the spectrum 



- ^=1 „ ^C T P gq (x)[l + (1 - xf - x 2 /9]Lne 2 p*(L) 2 / 1 

~ 2C F ° g \e 2 p 2 d (L), 



To calculate the N = 1 term in the OA one should replace in (11) l/(k 2 + mjy) 2 by 
2g5(k)/na 2 C < ACj' . In this case the N = 1 contribution vanishes in the massless limit 
[12, 15, 16], and all the contribution comes from the mass correction 

dP^ 4qLa s P gq (x)[l + (1 - x) 2 - x 2 /9]e 2 pj(L) ( 1 \ 

^c~f g VWL))- {) 

Thus, one sees that in the deep diffusion regime the gluon yield has anomalous mass 
dependence both for the realistic Hamiltonian and in the OA. 

In terms of the representation (3) the different mass dependence in the diffusion and 
infinite medium regimes is due to qualitatively different character of the p-dependence of 
the in-medium LCWF in these two situations. In the infinite medium regime ip(x, p, z), 
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like i/)(x, p), is a smooth exponentially decreasing function of p. In this case the probability 
of gluon emission decreases with quark mass due to a reduction of the dominating p scale 
in (3). Naively one could expect that the same mechanism should suppress the gluon 
emission from heavy quarks in the diffusion case. However, this is not true since in the 
diffusion regime the in-medium LCWF entering the integrand of (3) becomes oscillating in 
p. It is well seen from Fig. 1 where we plot the leading order radial in-medium LCWF for 
different values of the dimensionless parameter £ = L/Lf. One can see that the FS effects 
becomes small only at L/Lf > 10. The oscillations of the in-medium LCWF suppress 
the probability of gluon emission. This suppression is weaker for heavy quarks (since the 
ratio L/Lf is smaller). If the effect of the FS oscillations overshoots the mass suppression 
of the integrand in (3) which arises from the ordinary LCWF, the gluon spectrum rises 
with quark mass. Our calculations above demonstrate that namely this occurs in the deep 
diffusion regime. 

It is not surprising that the dead cone arguments of [7] fail in the diffusion regime since 
they are probabilistic in nature, while the oscillations of the in-medium LCWF leading 
to the anomalous mass dependence is a purely quantum effect. It is worth noting that 
the upper bound on the gluon momentum up to which the derivation of the dead cone 
suppression is supposed to be valid has been obtained in [7] from a crude estimate. In 
reality, as will be seen from our numerical results, the applicability region of the infinite 
medium approximation turns out to be considerably narrower, and in a broad kinematical 
range of gluon energy instead of the dead cone suppression we have an enhancement of 
the gluon yield. 

Note that, similarly to the gluon emission, the photon radiation from massive quarks 
is also enhanced in the diffusion regime. The N — 1 contribution in the massless limit has 
been calculated in [19] (similarly to the gluon spectrum it is oc L 2 ). The mass correction 

dP q N Zl _ e 2 q a em a 2 s C T C F nL 3 m 2 q x[l + (1 - x) 2 } 2 ( \_ 

dx ~ 16£ 2 (1 -x) 2 ° g \e 2 p 2 d (L) 

where now e 2 = m 2 x 2 . Due to large cross section of the charm production the photon 
radiation from c-quark may become an important mechanism of the photon production 
at LHC energies. Since the radiated photon has large momentum (x ~ 1) measurements 
of the photon tagged charm production would be especially interesting. 
4. The above analysis is very qualitative, and is not valid in the intermediate region 
L/Lf ~ 1. For an accurate evaluation of the mass dependence numerical calculations are 
necessary. We have performed computations (with any number of rescatterings) using the 
method suggested in [4]. It reduces calculation of the spectrum to solving the Schrodinger 
equation for the Hamiltonian (1) with a smooth boundary condition. As in [5] for light 
quark we use the quasiparticle mass m q = 0.3 and for gluon m g = 0.4 GeV obtained in 
[20] from the lattice data within the quasiparticle model. For the Debye mass we take 
rriD = \[2m g . For heavy quarks we take m c = 1.5 and = 4.5 GeV. We have evaluated 
the ratio of the gluon spectra for heavy and light quarks for plasma temperature T = 250 
MeV and a s = 0.4. In Fig. 2 we plot the results for E = 20, 50 and 100 GeV, and for the 
plasma thicknesses L = 2, 4 and 6 fm. The value L = 2 fm approximately corresponds 
to the typical opacity of the plasma produced at RHIC energies. The solid and dashed 
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curves show the results with and without the FS effects. One can see that at E > 50 
GeV the FS effects are very important at L ~ 2 — 4 fm. From the point of view of the jet 
quenching the soft region x < 0.5 is especially important. Fig. 2 shows that in this region 
for L ~ 2 fm and E ~ 100 GeV the FS effects increase considerably gluon emission from 
heavy quarks. The FS effects become small only at E ~ 20 GeV for L > 4 fm. Thus, 
Fig. 2 clearly demonstrates importance of the FS effects in the mass dependence of the 
gluon yield at L < 4 fm. 

For comparison with [7] we have also performed numerical calculations in the OA when 
the rescatterings can be characterized by the transport coefficient q. We take q = 0.3 
GeV 3 . At this value of q the gluon spectrum in the Bethe-Heitler limit x — > agrees with 
that calculated with accurate three-body cross section. In Fig. 3 we present the results 
for two sets of the parton masses. The solid curves correspond to the same masses as in 
Fig. 2, and the dotted ones are for massless gluon and light quark as in [7]. The dashed 
curves show the Dokshitzer-Kharzeev [7] dead cone suppression factor. Comparison of 
the dashed and dotted curves show that the results of the accurate calculations differ 
drastically from that predicted in the dead cone model. The results for massive and 
massless gluon in Fig. 3 show that the nonzero gluon mass (the Ter-Mikaelian effect) 
reduces the anomalous mass dependence. The anomalous mass dependence in Fig. 3 is 
stronger than in Fig. 2. It is quite natural since in the OA the leading N = 1 rescattering 
term simply vanishes for massless partons, contrary to the case of realistic Hamiltonian 
when due to the Coulomb effects the single scattering contributes even in the massless 
limit. 

From Fig. 3 one can conclude that the applicability region of the infinite medium 
approximation is considerably narrower than it was assumed in [7]. The authors of [7] 
assume that their model is valid for the gluon energy oo < ql? . This means that for L ~ 5 
fm the mass dependence can be described in terms of the dead cone suppression in the 
whole kinematical range of x at E < 100 GeV. Our results shown in Fig. 3 demonstrate 
that this is not true. Indeed, one can see that the FS effects come into play very early. 
As a result for x < 0.5, E ~ 100 GeV there is no the dead cone suppression at all. 

Our results indicate that there must not be a considerable difference between the jet 
quenching for c and light quarks at E ~ 20 — 50 GeV. This energy interval is interesting 
from the point of view of the nuclear modification factor Raa for the non-photonic elec- 
trons measured at RHIC [8]. If charm dominates in the region p? < 8 GeV studied in 
[8], then the observed suppression of the non-photonic electrons is not surprising. Even 
if at p T ~ 5 — 8 GeV about 50% of the electrons come from bottom quark [21, 22] it is 
hardly possible to speak of a serious contradiction of the theory and experiment. Indeed, 
in this scenario one can obtain the electron suppression about 0.3-0.5 while the data [8] 
give Raa ~ 0.1 — 0.55 with large error bars. More accurate comparison of the energy loss 
models with experiment will be possible after appearance of the data on the open charm 
production that are expected soon [23]. 

Note that numerical calculations within the path integral approach performed in [24] 
also indicated that in some parameter range the average energy loss grows with quark 
mass. The authors interpreted this fact as an unphysical effect related to inapplicability 
of the high-energy approximation. However, the situation with applicability of the LCPI 
formalism in the kinematical region dominating the anomalous mass dependence is in 
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fact quite good. Indeed, the LCPI formalism [9] is derived assuming that the longitudinal 
momenta are large compared to the transverse ones and the parton masses. Both these 
conditions are well satisfied in the diffusion kinematical region related to the anomalous 
mass dependence. In this regime the situation with the validity of the LCPI approach is 
even better than in the infinite medium regime (which gives the mass suppressed gluon 
yield) where the radiated gluons are softer. 

5. In summary we have shown that for a FS QGP the induced gluon radiation from 
heavy quarks becomes stronger than that for light quarks when the gluon formation 
length becomes comparable with or exceeds the size of the plasma. In this regime the 
gluon yield is dominated by the N — 1 rescattering. Physically the anomalous quark mass 
dependence is due to oscillations of the LCWF for the in-medium q — > gq transition. The 
dead cone model [7], which neglects the quantum FS effects, is not valid in this regime. 
The anomalous mass effect becomes more pronounced if one neglects the Coulomb effects 
and describes the parton rescatterings in terms of the transport coefficient. The neglect 
of the gluon quasiparticle mass also enhances the anomalous mass dependence. 

Similarly to the gluon emission, the FS effects enhance the photon radiation from 
heavy quarks. For this reason the photon radiation from c-quark may be an important 
mechanism of the hard photon production at LHC. 
Acknowledgements 

BGZ thanks the LAPTH for the kind hospitality during the time when part of this work 
was done. The work of BGZ is supported in part by the program SS-3472.2008.2 and the 
LEA Physique Theorique de la Matiere Condesee. 

References 

[1] P.M. Jacobs, M. van Leeuwen, Nucl. Phys. A774, 237 (2006) [arXiv:nucl-ex/0511013] 
and references therein. 

[2] R. Baier, D. Schiff, and B.G. Zakharov, Ann. Rev. Nucl. Part. 50, 37 (2000) 
[arXiv:hep-ph/0002198] and references therein. 

[3] A. Kovner and U.A. Wiedemann, arXiv:hep-ph/0304151 and references therein. 

[4] B.G. Zakharov, JETP Lett. 80, 617 (2004). 

[5] B.G. Zakharov, JETP Lett. 86, 444 (2007) [arXiv:0708.0816 [hep-ph]]. 

[6] B.G. Zakharov, JETP Lett. 88, 781 (2008) [arXiv:0811.0445 [hep-ph]]. 

[7] Y.L. Dokshitzer and D.E. Kharzeev, Phys. Lett B519, 199 (2001). 

[8] S.S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 96, 032301 (2006); J. 
Bielcik et al. [STAR Collaboration], Nucl. Phys. A774, 697 (2006). 

[9] B.G. Zakharov, JETP Lett. 63, 952 (1996); ibid 65, 615 (1997); 70, 176 (1999). 

[10] B.G. Zakharov, Phys. Atom. Nucl. 61, 838 (1998). 



7 



[11] B.G. Zakharov, Proceedings of the 33rd Rencontres de Moriond: QCD and High 
Energy Hadronic Interactions, Les Arcs, France, March 21-28, 1998, pp. 465-469; 
arXiv:hep-ph/9807396. 

[12] B.G. Zakharov, JETP Lett. 73, 49 (2001). 

[13] P. Aurenche and B.G. Zakharov, JETP Lett. 85, 149 (2007) [arXiv:hep-ph/0612343]. 

[14] N.N. Nikolaev, G. Piller, and B.G. Zakharov, J. Exp. Theor. Phys. 81, 851 (1995) 
[arXiv:hep-ph/9412344]; Z. Phys. A354, 99(1996) [arXiv:hep-ph/9511384]. 

[15] P. Aurenche, B.G. Zakharov, and H. Zaraket, JETP Lett. 87 605 (2008) 
[arXiv:0804.4282 [hep-ph]]. 

[16] P. Arnold, arXiv: 0903. 1081 [nucl-th]. 

[17] R. Baier, Y.L. Dokshitzer, A.H. Mueller, S. Peigne, and D. Schiff, Nucl. Phys. B483, 
291 (1997); ibid. B484, 265 (1997); R. Baier, Y.L. Dokshitzer, A.H. Mueller, and 
D. Schiff, Nucl. Phys. B531, 403 (1998). 

[18] B.G. Zakharov, JETP Lett. 80, 67 (2004). 

[19] B.G. Zakharov, JETP Lett. 80, 1 (2004). 

[20] P. Levai and U. Heinz, Phys. Rev. C57, 1879 (1998). 

[21] A.D. Frawley, T. Ullrich, and R. Vogt, Phys. Rep. 462, 125 (2008) [arXiv:0806.1013 
[nucl-ex]]. 

[22] A. Mischke, Phys. Lett. B671, 361 (2009) [arXiv:0807.1309 [hep-ph]]. 

[23] S.L. Baumgart [for STAR Collaboration] arXiv:0805.4228 [nucl-ex]; private commu- 
nication. 

[24] N. Armesto, C.A. Salgado, and U.A. Wiedemann, Phys. Rev. D69, 114003 (2004). 



8 



Figures 



Re(p ¥(p,Q) 



^> 1 1 1 1 1 


A ' A 


^~ 00 \ xT**^!^ \ / 

- — ^ =1 ° \ K 1 T 

-— ?=2 \ j\ \ 1 

-—5=i \ / \ 1 / 

- — ^ =0 - 5 \ / \ \ / 
— ?=04 V , H/ 




,V 1/V 11 



1 2 3 4 5 

P 



Figure 1: The leading order radial LCWF for in- medium q — > gq transition as a function of 
p in units of 1/e for different values of the dimensionless longitudinal coordinate £ = z/Lf. 
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Figure 2: The ratio of the gluon spectra for heavy and light quarks evaluated with the 
dipole cross section for the Debye-screened potential. The blue curves are for c-quark, 
the red curves are for 6-quark. The solid and dashed curves show the results with and 
without the FS effects. All the results are for m g = 0.4 GeV. 
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Figure 3: The ratio of the gluon spectra for heavy and light quarks evaluated in the OA. 
The blue curves are for c-quark, the red curves are for 6-quark. The solid curves show 
our results with FS effects for m g = 0.4 GeV, and the light quark mass m q = 0.3 GeV, 
and the dotted curves show our results for massless gluon and light quark as in [7]. The 
dashed curves show the Dokshitzer-Kharzeev dead cone suppression factor (9). 
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